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ABSTRACT 

We calculate the angular power spectrum of the galactic synchrotron radiation induced by the small scale 
fluctuations of the magnetic field and the cosmic ray electron density. Using the observed interstellar magnetic 
field spectrum, which is consistent with the Komolgorov turbulence model at the relevant scales, we find that 
Ci oc r^ '' . We estimate the cosmic ray electron density fluctuation spectrum with an injection-diffusion model, 
the shape of the angular power spectrum in this model depends on the correlations between the injection 
sources. For Poisson distribution of sources, C; oc /~^. We discuss the implications for the interpretation of 
cosmic microwave background (CMB) data and the impact on future 21 cm tomography experiments. 
Subject headings: cosmic microwave background — radiation mechanism: non-thermal — radio continuum: 
galaxies — ISM: magnetic fields — cosmic rays 



1. INTRODUCTION 

The radio sky at < IGHz is dominated by galactic 
synchrotron emission. It is believed to be produced by 
cosmic ray electrons propagating in the magnetic field of 
the Galaxy ( Ginzburg & Svrovatskii 1969). The galactic 
synchrotron emission is an important foreground for the 
cosmic microwave background (CMB) experiments dSmootl 
[T999). For the upcoming high redshift 21crn tomogra- 
phy experiments (Madau, Meiksin, & Rees'1997';'Tozzi et al.' 
2p00; Chen & Miralda-Escude 2004; Ciardi & Madau 2003; 
Gnedin & Shaver 2003; Furlanetto, Sokasian, & Hernauist 
.2004a; Iliev et al 2003; Loeb & Zaldarriaga 2004) , such as 
PAST jPen et aZ.l2004h . LOFAR' , and SKA^, it poses a major 
challenge. 

The global galactic synchrotron emission spectrum from 
0.3 MHz to 408 MHz can be fit with a two component 
disk model of the galaxy JK eshet. Waxman."& Loeb 200 4a,lj) . 
Analysis of the existing rad io surveys at 4 08 MHz , 1 .42 
GHz, and 2.326 GHz (Haslam ef a/.' 1982; Reich 1981 
|Reich& Reich 1986, 1988; Jonas. Baart. & Nicolson 199i 
shows that the synchrotron emission has a spectral index 
/? « 2.7 (Platania et al. 2003), which is in general agreement 
with the CMB result (Bennett et al. 2003). The real space dis- 
tribution of synchrotron emissivity over the galactic disk, tak- 
ing into account spiral arms, were derived from the 40 8 MHz 
whole sky map using refolding techniques ( Phillip s"ef al.\ 
01)81ab; Beuermann. Kanbach,& Bekhuiisen 1985). In the 
Fourier space, the angular power spectrum of the galactic syn- 
chrotron radiation can be modeled reasonably well as a sim- 
ple scaling relation: C; oc l~^ (Tegmark & Efstathiou 1996; 
[Tegmark.E isenstein.de Ohveira-Costal l2000t iGia rdino et al. 
UpOl, 2002), wifli 7y = 2.4-3 down to / ~ 900. Recendy, the 
WMAP team obtained a shallower spectrum of 77 ~ 2 down to 
/ ~ 200. These are extrapolated to higher / in recent studies 
of the 21 cm foregro und (Di Matteo. Ciardi. & Miniati 2004 
ISantos. Coorav. & Knox 2004). 

In the present study we investigate the galactic synchrotron 
emission from a different perspective. We calculate the syn- 
chrotron angular power spectrum induced by the variation of 
the magnetic field and fluctuations in the cosmic ray elec- 

' http://www.lofarorg 

^ http://www.skatelescope.org 



tron density. Our primary objective is to achieve a physi- 
cal understanding of the origin of the synchrotron emission 
anisotropy. In particular, we would like to ask which of 
these two mechanisms is responsible in producing the ob- 
served anisotropy? This will help us to assess the validity 
of the hypothesis adopted in the empirical analyses of CMB 
and 21 cm observation, e.g., will the power law form of the 
angular power spectrum hold down to small scales relevant 
for the pre-reionization 21 cm observation? At the same 
time, we may also gain useful knowledge about the galac- 
tic distribution of the cosmic ray electrons and the magnetic 
field. Similar physical modeling have been perf ormed for 
a number o f other fo regrounds, e.g. free-free ( Oh & Mac^ 
2003; Coorav & Furl anetto 2004.) and intergalactic shocks 
( Keshet et al. 2004a 1^. 

The fluctuation power spectr um of the galactic mag- 
netic field has been measured dMinter & Spangler 199^ 
iHan. Fe rriere. & Manchester'2004'). and is consistent with be- 
ing produced by a turbulent interstellar medium (ISM) de- 
scribed by the Komolgorov scaling model ( Komolgorov 194ll 
lElmegreen & Scaiol 12004ft . With this we can calculate di- 
rectly the anisotropy power spectrum induced by the magnetic 
field variation. The distribution of the cosmic ray electrons 
is less well known, but there is a broadly accepted picture 
of the cosmic electrons being produced in supernovae rem- 
nants (SNR), which then diffuses through the whole galaxy 
and be confined in a volume greater than the galactic disk 
(Ginzburg & Svrovatskii 1964). We can calculate the cosmic 
ray density distribution with this model. 

2. MODELS 

Let us consider the angular power spectrum obtained over 
a small patch of "blank" sky field at high galactic latitude. In 
the optically thin case, the observed radiation intensity is sim- 
ply an integral of the emissivity along the line of sight. At 
very low frequencies, synchrotron self-absorption and plasma 
absorption are important. For 70 MHz <v < 200 MHz, corre- 
sponding to the redshift range of 6-20 which is of interest to 
the study of the reionization process, these can be neglected. 



/(n,i^)= / dr e(h,r,iy), 



(1) 



where e is the volume emissivity. The emissivity may have 
both spatial and frequency variations. If the two are separable. 
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we may write 

e,{^) = g(um^), e,{^) = g{um^)- (2) 

where the tilde denotes Fourier transform, ^(x) = 
■ _£k_ 



/ W^^'^ '^''!^^^- The emissivity power spectrum is then 



(3) 



Let us consider the case v\=V2 = v. At small scale, using the 
Limber approximation ( Kaiser 1992; Zal darriaga et a/. 2004I) . 
the angular power spectrum is given by 



rT, . c'^\v) fdr I 



(4) 



If the emissivity has a power law spectrum, Ptpik) oc k^, then 
C/ oc The integral is up to some cutoff point, at which the 
emissivity drops to 0. 

We now consider the emissivity of synchrotron radiation. If 
the energy distribution of cosmic ray electrons f(E) at each 
point is approxiamated as a simple power law with f(E) = 
CE-P, then 
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where a,,„(p) = ^r[^]r[^]. We see that for power 

law distribution with f(E) oc E~'\ I^u) oc v"', and T{v) oc v^, 
where a = -(p- 1)/2, and (3 = a-2. A power law is actually a 
good approximation to the real distribution function of cosmic 
ray electrons. In this approximation, variation of emissivity 
can be induced by varying the magnetic field B, spectral index 
p, and cosmic ray electron density normalization C. Here we 
shall fix p and consider the variation of B and C. 

The galactic magnetic field is typically a few /iGs, for this 
magnetic field the emission at 70 -200 MHz is produced pri- 
marily by cosmic electrons with £ ~ 0.1 GeV- 10 GeV. The 
local cosmic ray electron density can be measured directly. 
There is some uncertainty in the normalization, and correc- 
tions have to be made for sola r modulation. A recent compi- 
lation of measurements yields dCasadei & Bindi2004h . in our 
units, C = 1.7 X 10"''cm"^GeV"', and p = 3.44, at the energy 
of a few GeV. The spectral index p = 3.4 yields a = -1.2, 
and /3 = -3.2. For comparison, the WMAP measurement 
(^ennett ef fl/. 2003) indicates that towards star forming re- 
gions, f3 = -2.5, and towards the halo, /3 = -3.0. Radiative 
loss of electron energy provides a natural explanation to the 
steeping of the electron spectrum away from star-forming re- 
gion. Below 3 GeV, the electron spectrum becomes flatter, 
probably because the primary energy loss mechanism changes 
from radiation to ionization. The locally measured electron 
density may not represent the average, nevertheless we will 
use it as a trial value. We will take p = 3, which gives a 
slightly better fit to the radio data than the steep /:> = 3.4 value. 
We also assume a disk scale height of 1 kpc (correspond- 
ing to the thick disk in the iKeshet et a/.ll2004al model), and 
a smooth magnetic field of 4 /xGs. With this set of parameters, 
the integrated sky brightness temperature at 408 MHz is 20 K, 
which reproduces the observed value at high galactic latitude 
(|Haslam ef fl/. 1982). 

Let us consider first the magnetic field with C fixed. Now 
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Fig. 1. — The magn etic field energy spec trum as extrapolated from 
IMinter & Soangleil 119961) and lHan ef oil 120041) . The dotted part is inter- 
polation. 



If there is a large smooth global magnetic field component 
which varies only on large scale and fluctuations around it are 
small, then we can write B = Bq + 6B(x), and 



B 



then 



P4,(k) ■■ 



Bqa 



C'B'o^^PsBik). 



(7) 
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The interstellar medium is turbulent ( lElmegreen & Scalol 
l2004h . Komolgoro v derived a scaUng relation for scale- 
invariant turbulence (iKomol ^orovl 1 94 ll) . with a power spec- 
trum of the form k~^^^^. The magnetic field fluctuation spec- 
trum can be determined by energy equipartition. The pre- 
dicted magnetic field fluctuation spectrum is confirmed on 
the scales of 0.01 pc - 100 pc by observation of the Farada y 
rotation of extragalactic sources (^ Minter & Spangle3ll996^ . 
At larger scales and on the galactic disk, the magnetic field 
spectrum is flatter, probably because at these scales the mo- 
tion is dominated by by two-dimensional structure (vortices). 
If we join the s mall scale and the large scale observations 
JHan et a/.l2004 . we obtain 

i 2.03 X 10-" erg cm'^ kpc, k^^-, > 1.57 x 10^ 
I 1.34 X 10-12 k^^Jleig cm-3 kpc, fe^pc-' < 6.28 

(9) 

This is plotted in Fig.[n(Note that our definition of k differs 
from Han_eLfl/. (2004) by a factor of 2tt). The energy is re- 
lated to the power spectrum by 

PB{k) = 2EB(k)/e. (10) 

Other scaling models have also been suggested. For exam- 
ple, some MHD turbulence have E ^ k~^/^ instead of k~^l^ 
(Kraichnan 1965). However, this would produce an angular 
power spectrum not very different from the Komolgorov one. 

We can then carry out the calculation, with the magnetic 
field spectrum given in Eq. (|9j and the locally measured elec- 
tron density. However, even though we keep C fixed to inves- 
tigate the variation induced by the magnetic field, in reality 
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Fig. 2. — The synchrotron power spectrum for v = 150MHz (upper curves) 
and 1/ = 20 GHz (lower curves). The solid line marks the anisotropy induced 
by the fluctuation of magnetic field, the dashed line for that induced by the 
fluctuation of cosmic ray electron density, the points are WMAP K band (23 
GHz) data. 



it must decrease as we move away from the galactic disk. 
We can approximate this effect by taking the power spec- 
trum Pti,(k) also as an explicit function of r, i.e. P^(k^r) ^ 
P{k)e~^''l''° , where ro is the scale height of the halo in which 
the cosmic ray electrons are confined. It turns out that this 
damping factor only changes the result by a small factor, be- 
cause the large r contribution is already suppressed by the 
P ^ (1/ r)"^/^ factor. For the small scale that we are mostly 
interested in, Komolgorov spectrum applies, Psik) ^ A;""/^, 
andC/ -r". 

The result of our calculation for v = 150 MHz (21 cm line 
at z ~ 9), and 23 GHz are plotted in Fig.|2]as solid lines, along 
with the WMAP K band data (centered at 23 GHz). Remark- 
ably, at small /, the anisotropy power amplitude is of the same 
order of magnitude as the WMAP data (At such small I, the 
Limber approximation may not be very good, nevertheless 
the true result will be of the same order of magnitude). This 
means that the magnetic field fluctuation may play a role in 
the formation of the observed synchrotron anisotropy. How- 
ever, it is clear that compared with the data C/ drops too fast 
as / increases. 

If the smooth magnetic field component is absent and the 
fluctuating field dominates, then 

^"^')(k) = {B"^')(k)c(p) 



B' 



(11) 



where c(p) = \ /^f dQ sin''-'^ 61 = r[(3 +/?) /2] /r[(2+p /2)] . 
In this case, the result is uncertain, because it depends on the 
higher order correlation {P-g{,M)ii{k)) . If we make the gaussian- 
like ansatz {BP^\k)) ~ ({B\k))yP^^'^/^, and Pgik) - r*?*, then 
the result is Q oc I'^-p+^'^'^bI^, For the Komolgorov spectrum 
TjB = -11 /3, C; ~ which is extremely steep and can be 
neglected entirely. However, this assumption maybe incor- 
rect. 

Now we consider the variation of the cosmic ray electron 
density. This is not known from observation. To make an es- 
timate of the fluctuation, we consider an injection-diffusion 
model (iGinzburg & Svrovatskiil[T964l) . in which the cosmic 



ray electrons are produced (injected) in point sources, then 
diffuse out, until eventually losing all their energy or escaping 
the confinement volume. The density of electrons at a point 
in space is then dependent on the distance to nearby sources. 
Neglecting the momentum space diffusion, the density nor- 
malization constant at position r, time t satisfies the diffusion 
equation 



5C(r,f) , 

—^=D^^C{v,t) + q{v,t)- 



C{r,t) 



(12) 



where q(r,t) is the cosmic ray electron injection rate at r, and 
D is the diffusion coefficient, which is constant bel ow 5 GeV, 
lO^\E/5GeV)°-^ at £ > 5GeV (Koba 



andP « 2 X 10^'*(£/5GeV)" *' at £ > 5GeV ( Kobavashi. 
120041) . The last term in the equation represents loss of elec- 
trons by radiation, ionization, etc., or by escaping the confine- 
ment volume, with r th e loss time scale. For radiative losses 
dCasadei & Bindil2004 . 

T w2.1 X 10^(£/TeV)-Vr. (13) 

We note that in the injection-diffusion model, the scale height 
of the cosmic ray halo is a few times of x VDt, with the above 
values we have \/Dt ~ 0.3kpc. Although very crude, our 
model is self-consistent. In Fourier space. 



■ + (Dk^ + -)C(k)-- 

T 



dt 

The steady state solution is 
C{k) = 

The power spectrum is then 
Pcik) 



■ qik) 



qik) 
Dk^+- 



(Dfc2+1)2' 



(14) 



(15) 



(16) 



Supernovae remnants are most likely the primary source for 
these cosmic ray electrons. The injection function is then 
q = NgKsN, where Ne is the number of cosmic ray electrons 
produced in one supernova, and ksn is the number of super- 
novae explosions per unit volume per unit time. If the distri- 
bution of supernovae is Poisson, with (ksn) = 1 /{VIsn), where 
V is the volume in which one supernova explode per average 
interval f^^v, then 

When ^ ^ 0, we have Pcik) = CqV, where Co is the average 
of C. From this we obtain 



V 



■■Co 



tsN 



(18) 



As a reality check, we take E ^ IGeV, a stellar disk with ra- 
dius of 15 kpc, and a scale height of 300 pc, and also assume 
that in this volume the average interval of supernovae explo- 
sion is 50 years, then we find Ne ^ 10"***, which requires an 
energy of 10"*^ erg, which is a small fraction of the total energy 
of a supernova, hence supernovae do have sufficient energy to 
generate these cosmic ray electrons. 
Using this relation we finally obtain 



Pr= V (19) 

Note that when written in this form, the result does not de- 
pend on tsN- The resulting C; is plotted as dashed curves in 
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Fig.|2] At large angle, the predicted C/ in this model is com- 
parable to the case of magnetic field induced fluctuation, and 
agrees with what is observed at the order of the magnitude 
level. However, at large I, Ci ~ again the spectrum is too 
steep compared with observation. Although the spectrum flat- 
tens at small k, this happened only at scales comparable to the 
disk scale height, thus affecting only large angle (/ ^ a few). 

In the above we have assumed a Poisson distribution of su- 
pernovae remnant. If they are correlated, with Pswik) oc Ak^, 
then Ci oc /'''"'*. Supernovae may well be correlated, as their 
rate should be proportional to the star formation rate, which 
in turn depends on the density. However, the Komolgorov 
model of turbulence suggests 7 < 0. The observation of floc- 
culent star light distribution in nearby spiral galaxies seem to 
confirm this expectation, which has a power law of Pa ^ 
at the 100 pc sca le (Elm egreen. E lmegreem. & Leitner 20031 
lElmegreen et aZ.l l2003h . If so, flien the correlation of SNR 
may not help us. However, more observations are needed to 
address this issue. 

3. CONCLUSION 

We have calculated angular power spectrum of the galactic 
synchrotron radiation induced by the variation of the magnetic 
field and cosmic ray electron density. We found that at low 
I, the amplitudes of the anisotropy power produced by both 
mechanisms are comparable to the observed value. This indi- 
cates that these physical mechanisms are relevant to the pro- 
duction of the observed foreground anisotropy. However, nei- 
ther of these two mechanisms c an produce the C/ ^ power 
spectrum observed by WMAP iBennett et al' 2003). and are 
also steeper than the older value of 7] 2.4-3. The magnetic 
field model induce power spectrum of the form C/ ^ , 
while for the electron density fluctuation it is C/ ~ /"^ if the 
spatial distribution of the SNR is Poisson. This may be reme- 
died if the spatial distribution of supernova remnants has the 
form PsNik) ^ k^- However, observations seem to indicate 
PsN(k)--k-K 



We made a number of simplifications in our calculation. We 
did not consider the detailed distribution of global magnetic 
field and cosmic ray sources in the Galaxy, nor do we con- 
sider accelerations outside SNR. It is unlikely that inclusion 
of any of these details would change our qualitative conclu- 
sion. Perhaps more important is our assumption of a universal 
energy spectrum of the electrons. In reality, the spectral in- 
dex varies from place to place, and this may induce additional 
anisotropy. To model this, we need to include momentum 
space diffusion. We plan to address these issues in subsequent 
studies. 

If the galactic synchrotron emission does have a steep an- 
gular power spectrum as suggested here, how could we rec- 
oncile this with the shallower power spectrum reported by the 
WMAP team? One intriguing possibility is that there may 
be another type of foreground, which has a frequency de- 
pendence similar to the synchrotron radiation at the relevant 
bands, and was mistaken as synchrotron radiation. Recently, 
several groups of researchers have suggested that some of 
the foreground identified as synchrotron by the WMAP team 
may actually be spinnin g dust (Lagache 2003; Finkbein^ 
2003; de Olivei ra-Costa e t al. 2004). This has a peak at 
10 GHz -20 GHz, which would not affect the 21 cm measure- 
ments. There might also be other foregrounds, e.g. of ex- 
tragalactic origin, which contributes to the radio survey. Al- 
ternatively, there may be other unknown mechanisms which 
produce the small scale anisotropy in the galactic synchrotron 
radiation. Further investigations are needed to identify what 
is responsible for producing the small scale anisotropy power 
For the two known physical mechanisms discussed here, the 
anisotropy power at small scales is much smaller than derived 
from simple extrapolation. 
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